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Abstract 



Wc have conducted deep JHKs imaging polarimetry of a ~8' x 8' area of the NGC 2071 star forming 
region. Our polarization data have revealed various infrared reflection nebulae (IRNe) associated with the 
central IR young star cluster NGC2071IR and identified their illuminating sources. There are at least 4 

■ IRNe in NGC2071IR and several additional IRNe are identified around nearby young stars in the same 
V/~j | field-of-view. Each illuminating source coincides with a known near-IR source except for IRS3, which 

■ is only a part of IRN2 and is illuminated by the radio source lc. Aperture polarimetry of each cluster 
source is used to detect unresolved circumstellar disk/outflow systems. Aperture polarimetry of the other 
point-like sources within the field is made in this region for the first time. The magnetic field structures 
(from ~1 pc down to ~0.1 pc) are derived using both aperture polarimetry of the point-like sources and 
imaging polarimetry of the shocked H2 emission that is seen as the dominant knotty nebulae in the Ks 

' band image; they are both of dichroic origin and the derived field directions are consistent with each other. 

The magnetic field direction projected on the sky is also consistent with that inferred from the 850 /im 
q ■ thermal continuum emission polarimetry of the central 0.2 pc region, but running roughly perpendicular 

I (~75°) to the direction of the large scale outflow. We argue that the field strength is too weak to align the 

c/3 . outflow in the large scale field direction via magnetic braking. 

Key words: circumstellar matter — infrared: stars — ISM: individual (NGC 2071) — polarization — 
stars: formation 
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J-j \ 1- Introduction resolution polarimetry with adaptive optics on current 

large telescopes (e.g., Perrin et al. 2004), and highly- 
SIRPOL is the polarimetry mode of the JHKs band si- polarized sources detected with the wide-field polarimetry 
multaneous infrared (IR) camera SIRIUS (Nagashima et could provide the best targets for future high resolution 
al. 1999; Nagayama et al. 2003) mounted on the 1.4 m telescopes at various wavelengths such as ALMA or 30- 
IRSF survey-dedicated telescope in South Africa, which m class optical-IR telescopes. In this paper, we present 
has been available since 2005 December. The instrument imaging polarimetry of the NGC 2071 star forming region 
is the first to provide deep and wide-field infrared polari- in Orion. 

metric images, and can in principle measure polarizations The NGC 2071 star forming region lies in the north- 
of all the 2MASS-detected sources within a field-of-view ern part of the Orion B molecular cloud (Lynds 1630) at 
of T.7 x 7.'7 in the JHKs bands simultaneously with 1% a distance of 390 pc (Anthony-Twarog 1982). A dense 
polarization accuracy. Such a capability is unique and has molecular and dust cloud "core" is observed (e.g., Takano 
vast and unique applications in many fields in astronomy, et al. 1986; Dent et al. 1989). There is a cluster of IR 
including efficient detection of circumstellar material that sources in this region, NGC 2071IR, and at least 9 near- 
produces scattering (see Weintraub, Goodman, & Akeson IR and mid-IR sources have been detected in the core 
2000; Tamura & Fukagawa 2005 for a review), and delin- region (Persson et al. 1981; Lada et al. 1991; Walthcr 
eating magnetic field structures in molecular clouds (e.g., et al. 1993). A red near-IR source (SSV41) originally 
Tamura et al. 1987) and in galaxies (e.g., Jones 2005), found with a single-detector raster scan of a 31" beam 
both through dichroic polarization. by Strom, Strom, & Vrba (1976) seems to correspond to 

Wide-field polarimetry is complimentary to the high some of these collectively. Snell & Bally (1986) detected 
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three radio continuum sources, some of which (NGC 2071 
la and lb) are likely to have counterparts to the near-IR 
sources. The total infrared luminosity of the core region 
is estimated to be ~520 L Q (Butner et al. 1990), which 
suggests formation of intermediate-mass stars. 

An energetic bipolar CO outflow is centered near 
NGC2071IR, lying in the NE-SW direction (see Moriarty- 
Schieven, Hughes, & Snell 1989 and references therein). 
Molecular hydrogen emission is also detected (Eisloffcl 
2000; Aspin, Sandell, & Walther 1992 and references 
therein), which is well correlated with the CO outflow 
morphology. Dense gases associated with either outflows 
or disks are observed in CS, SO, SiO, HCO+, and NH 3 
(Kitamura et al. 1992; Chernin & Masson 1993; Girart 
et al. 1999; Zhou, Evans, & Mundy 1990). Several water 
masers and an OH maser are detected (e.g., Sandell et al. 
1985). Water maser emission mapping with VLB A (Seth 
et al. 2002) suggests that IRS1 and IRS3 are the exciting 
sources and both sources have compact (—10 AU) disks. 

The region has been one of the first targets of IR po- 
larimetry. Aperture polarimetry of IRS1 and its mapping 
were conducted by Heckert & Zeilik (1984). They de- 
tected a relatively large (~7% in the K band at 0=142°) 
polarization perpendicular to the outflow direction with a 
large 32" beam. Hodapp (1984) and Sato et al. (1985) 
suggested the presence of IR reflection nebulae (IRNe) 
extending roughly in the direction of the large scale out- 
flow. Therefore, the aperture polarization with a large 
beam can be explained with the dust scattering in the 
IRN. More recently, Walther et al. (1993) presented high 
resolution (0."62 pixel -1 ) K band imaging polarimetry of 
a -65" x 40" region of NGC2071IR, including IRS1 to 
IRS8a and b, as well as JHK & NBL band imaging of 
a ~40" x 35" region and K band imaging of a 5' x 9' 
region. However, no wide-field, deep imaging polarimetry 
or JH band polarimetry has been conducted yet. 

Information regarding magnetic fields in the NGC 2071 
region was obtained from 850 /im polarimetry (Matthews, 
Fiegc, & Moriarty-Schicven 2002), which suggests the 
presence of a weakly aligned magnetic field over —2 pc 
scale whose direction is roughly perpendicular to the out- 
flow direction. A very large scale (15° ~ 100 pc) mag- 
netic field structure including the Orion molecular cloud 
complex was studied with optical interstellar polarimetry 
(Appenzcllcr 1974), but their field star polarimetry data 
are too sparse to discuss the magnetic field structure in 
the NGC 2071 region (< 1 pc scale). 

2. Observations 

Observations of the NGC 2071 region were made on 
2006 February 23. The 1.25 (J band), 1.63 (H band), 
and 2.14 (Ks band) /im imaging polarimetry data were 
obtained simultaneously with the SIRIUS camera and its 
polarimcter on the 1.4-m IRSF telescope in South Africa. 
The SIRIUS camera can provide JHKs images simulta- 
neously, with a field of view of T.7 x 7.'7 and a pixel scale 
of 0."45 pixel™ 1 . The camera is equipped with three 1024 
x 1024 HgCdTe (HAWAII) arrays, three broad-band fil- 



ters, and two beam splitters for the simultaneous imaging. 
The polarimcter is composed of an achromatic (0.5 - 2.5 
(im) wavc-platc rotator unit and a high-extinction-ratio 
polarizer, both of which arc located upstream of the cam- 
era and at room temperature. See Kandori et al. (2006a) 
for the details of the polarimcter. 

The total integration time was 900 s per wave-plate an- 
gle. We performed 10 s exposures at 4 wave-plate angles 
(0°, 45°, 22°. 5, and 67°. 5) at 10 dithered positions (1 set). 
The same observation sets were repeated toward the ob- 
ject (9 times) and sky coordinates. The sky condition was 
excellent and the seeing size during the observations was 
0"9 (2 pixels) in the J band. 

After image calibrations in the standard manner using 
IRAF (dark subtraction, flat-fielding with twilight-flats, 
bad-pixel substitution, and sky subtraction), the Stokes 
parameters (I, Q, U), the degree of polarization P, and the 
polarization position angle 9 were calculated as follows: 

Q = Iq — I45, U = 1-22,5 — ^67.5, 

I = (Jo + h 5 + 122.5 + h?. 5 )/2, 

P = ^Q 2 + U 2 /I, and 9 = (1/2) arctan([//Q) . 

The accuracy of polarization position angle (59) was es- 
timated to be better than 3° based on the measurements 
during the commisionning run. The polarization efficien- 
cies at JHKs are also higher than 96%; thus no corrcetion 
was made on the polarization degree (P). 

3. Results and Discussion 

3.1. Intensity and polarization properties in the overall 
NGC 2011 region 

Figures 1 (a) and (b) show the JHKs color-composite 
I image and the JHKs color-composite PI image in two 
different scales, respectively. Figure 1 (c) summarizes the 
identification of notable features discussed in this section. 

To the southern part of the I image, there is diffuse 
near-IR emission associated with the optical reflection 
nebula NGC 2071. The illuminating source of the neb- 
ula, HDE 290861, is just outside of the image. 

The large, asymmetric bipolar nebula is centered in the 
field (~4' north of the optical nebula), in which the most 
striking feature is reddish knotty emission in the Ks band. 
From the similarlity between our Ks band image and 
previous H 2 emission images (Eisloffel 2000; Aspin et al. 
1992; Garden, Russell, & Burton 1990), we consider that 
the emission in our Ks band image comes mostly from 
the H2 emission. The Ks band emission continuously ex- 
tends from near IRS1 to the north-cast direction. There 
is a gap of the Ks band emission in the south-west, and 
knotty H2 emission again appears —3' to the south-west 
(NGC 2071SW; Walter et al. 1993). 

In contrast, JH emission is more smoothly distributed 
and forms a U-shaped nebula centered near IRS1 extend- 
ing to the north-east. It traces some arc structure in the 
southern part of the NGC 2071SW region. There is also 
a gap of the JH emission to the south-west. 
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Although both of the Ks band emission and the JH 
band emission show the asymmetric bipolar structure, 
their detailed distributions are different, as seen from the 
color difference in the / image. If we overlap these two 
emission components, the bipolar structure is more clearly 
recognized (as in figure la), and their distribution is very 
well correlated with that of the molecular bipolar out- 
flows such as the HCO+ outflow (Girart et al. 1999). Two 
near-IR emission peaks in the north-east coincide with the 
blueshifted HCO + emission peaks (Bl and B2) and those 
in the south-west with the redshifted peaks Rl and R2. 
Therefore, it is natural to consider that the molecular out- 
flows and the near-IR emission or reflection nebulae in the 
core region are related to the large scale outflow activity. 
The high spatial resolution Fabry-Perot H2 observations 
of Aspin et al. (1992) indicate that IRS1 is the likely 
source for the large scale outflow (see also §3.2). 

The appearance of these nebulae recognized in the I 
images is drastically different in the PI image. First, the 
diffuse nebula associated with the optical reflection nebula 
NGC2071 disappears, suggesting that their polarization is 
low at near-IR wavelengths. This is due to the contribu- 
tion of near-IR continuum from small dust grains (Sellgren 
1984). Second, the Ks band nebulae associated with the 
bipolar outflow cannot be seen, suggesting that their po- 
larization is also low. This is apparently due to the con- 
tribution of H2 emission (see §3.3). Third and in contrast, 
several highly-polarized nebulae become distinct. There 
are at least 8 such nebulae in the / image. They are most 
likely reflection nebulae, thus referred as IRN(e) in figures 
1. Some of them (IRN4, 7) are not distinct in the / image. 

The illuminating sources for IRN1-4 in the core region 
are discussed in detail below. Those for IRN5, IRN6, and 
IRN9 are SSV37 (Strom et al. 1976), IRS18 (Walther et 
al. 1993), and SSV40 (Strom et al. 1976), respectively. 
Those for IRN7 and IRN8 are new point-like sources seen 
in our intensity images. 

Nine IR sources (IRSl-8a, and 8b) in the core region 
observed by Walther et al. (1993) have been confirmed 
in our three-color composite image of the central cluster 
NGC2071IR although some of them are simply a part of 
near-IR nebula and not distinct in our images (see figure 
2 for their identification). Their relative positions, and 
aperture photometry and polarimetry are summarized in 
Table 1. Note that the errors in the position angles (9) 
include the uncertainty in the calibrations (see §2). The 
positional offset in the table is relative to the 2MASS posi- 
tion of IRS1 (05 /l 47 m 4. s 77, 00°21'42."8; J2000) with an ac- 
curacy of ~0"1. Photometric calibration was made with 
using the 2MASS-PSC. 

There are many other point-like sources in our intensity 
images which have not been reported before. We did not 
identify or list each source in this paper, but the calibrated 
FITS images are available upon request. 

3.2. Polarization images of the central NGC2071IR re- 
gion 

Figure 2 shows the KsHJ polarization vector maps of 
the central 61" x 45" region. Identifications with known 



IRS are also indicated. 

There are several IRNe associated with the central star 
cluster NGC2071IR. As has been suggested by Walther 
et al. (1993), prominent IRNe in the Ks band are those 
around IRS1 and IRS8 with a centro-symmetric vector 
pattern. We also note another prominent pattern around 
IRS4. Our Ks band polarization vector map is mostly 
consistent with the K band polarization vector map of 
Walther et al. (1993). However, there are both similarities 
and differences between those two vector images in details, 
which we compare below. A summary of these discussions 
is given at the end of this subsection. 

3.2.1. IRS1 

The vector pattern around IRS1 is consistent between 
the two data. IRS1 itself is highly polarized (26% at 
Ks) and the polarization vectors show a gradual direc- 
tion change over this source (like the character "S" from 
north to south). A large part of the IRN associated with 
NGC2071IR (whose extent is more than 30") is illumi- 
nated by IRS1, thus we refer to it as IRN1. This corre- 
sponds to the white nebula in the PI image (figure lb). 
The levels of polarization degree of IRN1 range from 25% 
to 45% at Ks. IRS1 positionally coincides with the radio 
continuum source NGC 2071 la (Snell & Bally 1986) and 
is considered to be self-luminous. 

The brightest part of IRN1 is situated just to the east of 
IRS1, whose direction is rather oblique to the large scale 
outflow direction. However, the axis (79°) of the compact 
(~10 AU size) disk suggested from the distribution of H 2 
masers around IRS1 (Seth et al. 2002) is consistent with 
the main extension of IRN1 relative to IRS1. The position 
angle of our software aperture polarimetry (Table 1) of 
IRS1 (9 = 177° in the Ks band) is consistent with the 
orientation of the compact disk ((177° - 90°) = 87° w 
79°). The compact disk around IRS1 must be tilted so 
that the west side is nearer to us if we assume the typical 
morphological relationship between the outflows and IRNe 
(e.g., Tamura et al. 1991). 

Our Ks band magnitude of IRS1 is 0.8 mag fainter than 
that in Walther et al. (1993), although JHKs magnitudes 
for this source and other sources are roughly consistent 
with each other (< 0.5 mag). Note, however, this is not 
true for IRS5 and several very faint magnitudes (> 18 
mag). 

3.2.2. IRS2 

IRS2 is not significantly polarized in Walther et al. 
(1993), but in our data it is polarized (7-17%) at a level of 
the polarization of nearby nebulosity, especially at shorter 
wavelengths (see Table 1). IRS2 positionally coincides 
with the radio continuum source NGC 2071 lb (Snell & 
Bally 1986) and is considered to be self-luminous. 

There is a local intensity peak and a slight change of 
vector pattern to the east of IRS2, but it does not show 
a centro-symmetric pattern like IRS1 and IRS8. Thus 
the relationship of this local peak and IRS2 in not clear. 
Higher spatial resolution polarimetry is necessary to reveal 
these structures. Note however that a local H2 emission 
peak is associated with IRS2 (Aspin et al. 1992). We con- 
sider that IRS2 is a YSO without extensive circumstellar 
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material. 
3.2.3. IRS 3 

IRS3, a red source just north of IRS1, itself is highly po- 
larized (P = 19% in the Ks band) and the vector pattern 
around IRS3 is partly centro-symmteric, whose perpen- 
dicular lines apparently point toward the south-west. 

This position is in fact coincident with a radio source 
NGC 2071 lc (Snell & Bally 1986). Because of the limited 
spatial resolution in both radio and NIR observations, it 
is difficult to make absolute astrometry. However, if we 
assume IRS2, which is point-like and shows no centro- 
symmetric polarization pattern, coincides with the radio 
source lb, then IRS1 coincides with la, but IRS3 does not 
coincide with lc. Their positions are plotted in figure 2. 

Therefore, we suggest that IRS3 itself is not self- 
luminous but is illuminated by the radio source lc, which 
is situated just south-west of IRS3 and invisible even at 
near-IR wavelengths. Such suggestions of invisible illumi- 
nating sources are also made in other star forming regions 
(e.g., GGD27 NIRS2; Tamura ct al. 1991) 

The nebula associated with IRS3 is thus an indepen- 
dent, compact IRN besides the large IRN associated with 
IRS1. We refer to it as IRN2 (the red IRN in the PI 
image). The levels of polarization degree of IRN2 range 
from 11% to 33% at Ks. Their northern part, however, 
has a different pattern; this part is illuminated probably 
by IRS1. 

Seth et al. (2002) suggested the presence of a nearly 
edge-on circumstellar disk around IRS3 whose axis is in 
the direction of 35°. The relative position between the 
radio continuum source and IRS3 is consistent with the 
idea that IRS3 is a part of IRN illuminated by the invis- 
ible radio continuum source. If this is the case, the disk 
around IRS3 is tilted so that the south-west side is nearer 
to us. Similar to the case of IRS1, the position angle of 
the aperture polarization of IRS3 (0 = 124° in the Ks 
band) is consistent with the orientation of the compact 
disk ((124° - 90°) = 34° w 35°). 

Eisloffel (2000) suggested that IRS3 is the exciting 
source of the H2 emission (the flow IIA and IIB), extend- 
ing roughly east and west. However, we consider that 
their flow IIA is just a part of the outflow associated with 
IRS1 and the flow IIB is with IRS8. IRN1 and IRN4 
(see §3.2.8) are their counterpart IRNe, respectively. In 
fact, their narrow-band H2 emission images might suffer 
from the continuum emission; our imaging polarimetry of 
these components demonstrates that IRS3 (or its illumi- 
nating radio source lc; see above) is illuminating none of 
these. Aspin et al. (1992) could not see any convincing ev- 
idence for activity centered on IRS3 in their Fabry-Pcrot 
H2 imaging data. 
3.24. 

IRS4 itself is not highly polarized as found by Walther 
et al. (1993) and ours, but our data suggest a clear asso- 
ciation of a compact reflection nebula with a centrosym- 
metric pattern around IRS4, extending to the south-east. 
We refer to it as IRN3. 

Walther et al. (1993) suggested IRS4 to be a back- 
ground star. However, we consider this to be a member of 



the NGC2071IR cluster because of the association of the 
local nebula IRN3. IRS4 is also associated with faint H2 
emission (Aspin et al. 1992). 

3.2.5. IRS 5 

IRS5 seems to be a local peak of the IRN1 associated 
with IRS1. Walther et al. (1993) suggested IRS5 as well 
as IRS5a and IRS8a to be shocked H2 peaks. Our Ks 
and H band magnitudes for IRS5 are more than 0.5 mag 
fainter than those in Walther et al. (1993). 

3.2.6. IRS6 

We found that IRS6 is a binary with a separation of 2" 
and is not highly polarized. 

3.2.7. IRS7 

IRS7 is a rather isolated red-color source, and it is rel- 
atively highly polarized (^10%). It was suggested to be 
an independent, local outflow source by Walther et al. 
(1993). 

Our polarization vector map shows no clear ccntro- 
symmetric pattern around IRS 7. However, its position 
angle of polarization is completely different from that of 
nearby point-like sources (see §3.4), which is suggestive 
of the circumstellar structure such as disk or envelope 
(e.g., Tamura et al. 2006). Note also that CO band emis- 
sion is detected in IRS 7, inferring that this source is a 
young stellar object (YSO; Walther, Gcballc, & Robson 
1991). Therefore, we suggest that IRS7 is a YSO with a 
disk/outflow system. 

3.2.8. IRS 8 

IRS8 is a faint red source to the sout-west of IRS1 sug- 
gested to be a YSO (Walther et al. 1993). Associated 
is a near-IR nebula which shows a clear centro-symmetric 
pattern, thus we believe IRS8 to be a self-luminous illumi- 
nating source. We refer to the bipolar nebula around IRS8 
extending both to the north-west and to the south-east as 
IRN4. 

3.2.9. EAST and WEST 

The nebula "EAST" in Walther et al. (1993) is a part of 
the large IRN1 illuminated by IRS1. The nebula "WEST" 
in Walther et al. (1993) is a part of the IRN4 illuminated 
by IRS8. The small counterpart of this IRN is also seen 
just to the south-east of IRS8. 

3.2.10. Summary 

We have identified four infrared reflection nebulae 
(IRN1-4) and their illuminating sources in the core re- 
gion (the NGC2071IR cluster). IRN1 is the brightest part 
of the IR nebula associated with the large scale molecu- 
lar outflow. Ncar-IR source IRS3 is simply a part of the 
nebula IRN2. IRN1, IRN2, IRN3, and IRN4 are illumi- 
nated by IRS1 (radio continuum source la), radio contin- 
uum source lc, IRS 4, and IRS8, respectively. The vec- 
tor pattern is centro-symmetric around IRS1, 4, and 8, 
and partly centro-symmetric around lc. No distinct local 
polarization pattern is seen around IRS5 or IRS6. IRS7 
has an intrinsic polarization, suggestive of an unresolved 
disk/outflow system. 

All the features discussed above are most clearly seen 
in the Ks band but are very similar to those seen in our 
H band vector image. Only a part of these IRNe is seen 
in our J band image because of extinction in the cloud. 
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3.3. Hi emission polarization and magnetic field struc- 
tures 

Another interesting feature of the nebula polarimetry is 
seen in the Ks band, in which the knotty H 2 emission is 
significantly polarized and their vectors are relatively well 
aligned with each other. Note that our Ks band / image 
is deep enough that most of the H 2 emission features de- 
tected by Eisloffel (2000) with narrow band imaging are 
recognized in our Ks band image (red knots in figure 1). 

Since the shocked H2 emission itself is not polarized in 
situ, its polarizations are likely due to dichroic absorption 
by aligned grains between the H 2 emission region and the 
observer, thus tracing the magnetic field structure of the 
NGC 2071 region. Such H 2 emission polarimetry has been 
conducted for only two other regions so far (Orion, Hough 
et al. 1986; Burton et al. 1991; Chrysostomou et al. 1994; 
DR21, Itoh et al. 1999). 

Figures 3 (a) and (b) shows the polarization vector map 
and the histogram of the position angle of the polarization 
for the brightest H 2 knots IA, respectively (see Eisloffel 
2000). Polarization vectors are shown only when their 
signal-to-noise ratio is better than 4. Therefore, the po- 
sition angle error of each vector is better than 8°. The 
dominant direction is PA ~ 115° (FWHM ~ 40°), which 
is roughly perpendicular to the outflow direction. 

Levels of apparent polarization degrees are relatively 
higher (> 10%), but this is probably due to our low signal- 
to-noise ratio, especially for the faint knots. The brightest 
knots have polarizations of a few percent, consistent with 
the level of dichroic polarization in this region discussed 
below. 

Figures 3 (c) and (d) shows the similar images for the 
H 2 knots V. Polarization vectors are shown only when 
their signal-to-noise ratio is better than 2. Because of the 
weaker surface brightness, the number of available pixels 
for useful polarizations are limited, but we can still see 
an alignment tendency. Interestingly, the dominant direc- 
tion is PA ~ 150°, different by —40° from the direction 
in the knots IA, but with a large dispersion (FWHM ~ 
90°). Therefore, the magnetic fields might be slightly bent 
near these knots. Deeper H 2 polarimetry is necessary to 
confirm this result. 

The southern H 2 knots (IB) are contaminated by re- 
flection or ionized nebulosity that is seen as an arc to the 
south. We believe this part to be illuminated by the cen- 
tral star of NGC2071 (HDE 290861) rather than by IRS1 
because only the edge facing to HDE 290861 is bright. 
Unfortunately, the signal to noise ratio in polarization is 
too small and two competing mechanisms make the po- 
larization pattern too complex to discuss either magnetic 
field structures or reflection structures. However, it is 
noteworthy that some of the H 2 knots IB also show the 
polarizations (not shown here) consistent with those in 
the knots IA. Thus, the magnetic field direction seems to 
be roughly the same in the northern and southern outflow 
regions. 



3.4- Aperture polarimetry and magnetic field structures 

No interstellar polarimetry has been made to decide 
the magnetic field structures in this ~1 pc region so far. 
Therefore, we measured software aperture polarizations 
of point-like sources detected in the field-of-view. We re- 
jected the sources whose photometric errors are greater 
than 0.1 mag or whose polarization is greater than 10%. 
Note that the position angle error of each vector is bet- 
ter than 7°. Such a large polarization is unlikely to be 
dichroic origin in this region (see §3.6). Also rejected are 
the bright saturated sources. 

A total of 53, 79, and 73 sources are measured in the 
J, H, and Ks bands, respectively. The aperture radius 
was 3 pixels. Aperture polarization vector maps (Ks) of 
point-like sources superposed on the / images are shown 
in figure 4. The distribution of position angles at JHKs 
is shown in figure 5. 

Although the dispersion is relatively large at shorter 
wavelengths, there is a distinct trend of their position an- 
gles; they are well aligned with each other and considered 
to be of dichroic origin (e.g., Tamura et al. 1987), which 
suggests a large scale magnetic field over the observed area 
(on ~1 pc scale). The dominant field direction is PA ~ 
120° with some deviation in the north-east region of the 
observed field. Note that the polarization of the central 
source of IRN6 (sec figure lb) has a clear intrinsic polar- 
ization, which is almost perpendicular to the polarization 
of other sources. 

There is a slight shift of the peak angle from Ks band to 
J band by ~10°. The dispersion is larger at shorter wave- 
lengths; the FWHMs of the histogram are ~20°, -30°, 
and ~50° in the Ks, H , and J bands, respectively. The 
origin of this wavelength dependence of position angles is 
unknown at present. 

This direction of the aperture polarization is consistent 
with that derived from the H 2 polarimetry (PA=115°). 
Both of these are due to dichroic absorption in the NGC 
2071 region and give the magnetic field direction projected 
on the sky. The only difference is that the aperture po- 
larimetry (of background stars) traces the entire cloud 
along the line-of-sight, while the H 2 polarimetry traces 
the front-half of the cloud. Since these directions are 
consistent with each other, both the aperture polarime- 
try of field stars and the H 2 emission polarimetry suggest 
that the magnetic field runs at PA=120° over the out- 
flow region, which is roughly perpendicular to the large 
scale outflow direction (~45°). This direction is some- 
what different from the very large scale field in the Orion 
region (Appenzellcr 1974) but rather consistent with the 
direction near the OMC-1 region (~120°; Kusakabe et al. 
2006). 

3.5. Alignment of outflow/disk systems with magnetic 
fields 

Alignment of an outflow with the magnetic field is dis- 
cussed in Matsumoto & Tomisaka (2004). They demon- 
strated a misaligned outflow, which is ejected in the direc- 
tion of the magnetic field in the 10 — 100 AU scale, while 
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its direction is inconsistent with the magnetic filed in the 
cloud scale. The alignment is controlled by the magnetic 
braking during the protostcllar collapse. The magnetic 
braking time is estimated as (see eq. [29] in Matsumoto & 
Tomisaka 2004), 



1/3 



V a ' 



(1) 



where Xj = c s (n /Gp) 1 / 2 and v a — B/^-np) 1 / 2 denote the 
Jeans length and the Alfven speed, respectively. The pro- 
tostcllar collapse is characterized by the free fall time, 
defined as, 
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The ratio between above two time scales is expressed as a 
ratio between the sound speed and the the AlfVen speed, 



h 
t s 



1/3 
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Matthews et al. (2002) estimated the magnetic field 
strength by examined their submillimeter data using the 
method of Chandrasekhar & Fermi (1953), and obtained 
the field strength of B 56^t G, assuming the rms velocity 
of Sv = 0.5kms _1 and the density of p = 3.8 x 10 -20 gcm~ 3 
(the number density of n = 10 4 cm -3 ). For the sound speed 
c s in equation (3), we should adopt the value of the rms 
velocity Sv, which is the so-called effective sound speed, 
and it includes the effects of the internal gas motion such 
as turbulence. Applying these values, equation (3) yields 
tb/tg = 2.4. This ratio indicates that the magnetic brak- 
ing does not have enough time to align the magnetic field 
of the protostar scale with respective to the cloud scale 
magnetic field in the period of the protostcllar collapse. 
In other words, the magnetic field is too weak to align the 
outflow with the magnetic field on the large scale. 

Although the ratio, tb/tg, estimated above is a some- 
what marginal value, it seems to be under-estimated. This 
is because the method of Chandrasekhar & Fermi (1953) 
has a possibility of over-estimation of the magnetic field 
strength as pointed out by Ostriker et al. (2001), Padoan 
et al. (2001), and Matsumoto et al. (2006). If the mag- 
netic field strength is weaker than 56/i G, the ratio, tb/tg, 
has a larger value. 

This is in contrast to the OMC-1 region and the 
NGC7538 region where the new born massive YSOs bend 
the local magnetic field near YSOs but the outflows tend 
to align in the outer (~0.5 pc) regions (Momosc et al. 
2001). 

The weakness of the magnetic field is also exhibited in 
the position angles of the aperture polarization of YSOs, 
since they are indicators of the local disk/outflow systems 
(see §3.2). Their position angles represent the position 
angle of disk, perpendicular to the outflow direction. The 
aperture polarizations of IRS1, IRS3, IRS7, and IRS8 have 
the position angles which arc different with each other 
(see Table 1). This large dispersion in the position angle 



also indicates that the magnetic braking is ineffective in 
aligning of the axes of YSOs in this region. 

However, it should also be noted that in clustered star 
forming regions such as the NGC 2071 region, the direc- 
tion of the disk/outflow systems could be affected by the 
other YSOs, which might result in a "random" distribu- 
tion of the disk/outflow direction. A systematic wide-field 
polarization survey toward various clustered star forming 
regions will be required for understanding the relationship 
between the YSOs' disk/outflow alignment with the local 
magnetic field. 

3. 6. Near-IR vs. submillimeter polarimetry 

It is interesting to compare the magnetic field direction 
derived from dichroic polarization at near-IR wavelengths 
and that from thermal emission polarization at submil- 
limeter wavelengths because Goodman ct al. (1995) 
claimed that the former was not always a good tracer of 
the magnetic fields within the dense clouds. They sug- 
gested a possible inefficiency of grain alignment in a dense 
region. 

Figure 6 shows a comparison of the PA histogram in the 
H band and that at the wavelength of 850 (im (Matthews 
et al. 2002). The 850 ^m emission traces the mag- 
netic field structure associated with the dense cloud of 
the NGC2071IR region (covering about 2' = 0.23 pc in 
diameter), while the near-IR data cover a ~1 pc scale re- 
gion. The general trend of the two directions shows a good 
agreement with each other, although their spatial scale is 
somewhat different. Therefore, our results demonstrate 
that both near-IR dichroic polarization and submillimeter 
emission polarization trace the magnetic field structures 
of ~0.1-1 pc scale associated with the NGC2071IR region. 

The agreement of the magnetic field structures seen at 
near-IR and far-IR/submillimctcr wavelengths is also re- 
vealed in the OMC-1 region (Houde et al. 2004; Kusakabe 
et al. 2006), the Ophiuchius dark cloud (Tamura et al. 
1996), and the NGC 2024 region (Kandori et al. 2006b) 

3. 7. Polarization efficiencies 

Figure 7 shows the degrees of polarization in the H 
band vs. the H — Ks color. There is a weak correlation 
as is expected from the dichroic polarization (§3.4). The 
average slope P{H)/{H — Ks) is 6.0±4.8, which is consis- 
tent with that (~5.6) in the M42 region (Tamura et al. 
2006; Kusakabe et al. 2006). Therefore, the polarization 
efficiencies have a similar value in the two distinct regions 
in the Orion cloud complex. Detailed comparison of the 
polarization efficiency with other star forming regions will 
be presented in a future paper. 

4. Conclusion 

We have conducted deep JHKs imaging polarimetry 
of a 7'7 x 7'7 area of the NGC 2071 star forming region. 
Main conclusions are summarized as follows. 

1. Our polarization data have revealed various infrared 
reflection nebulae (IRNe) associated with the cen- 
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tral IR young star cluster NGC2071IR and identi- 
fied their illuminating sources. There are at least 4 
IRNe in NGC2071IR. IRN1 is the brightest part of 
the IR nebula associated with the large scale molec- 
ular outflow. IRN1, IRN2, IRN3, and IRN4 are 
illuminated by near-IR source IRS1 (radio contin- 
uum source la), radio continuum source lc, near- 
IR source IRS 4, and near-IR source IRS8, respec- 
tively. We have suggested that IRS3 is simply a 
part of IRN2 illuminated by a (near-IR invisible) 
radio source NGC 2071 lc. Four additional IRNe 
(IRN5-8) arc identified around nearby young stars. 
The optical reflection nebula, NGC 2071 itself, is not 
distinct in the near-IR polarized intensity images. 

2. Aperture polarimetry of IRS1 and IRS3 shows that 
their polarization position angles are both consistent 
with the orientations of the compact (~10 AU) disks 
traced by water maser emission. 

3. Aperture polarimetry of the point-like sources 
within the field in the JHKs bands is made in this 
region for the first time. The magnetic fields (from 
~0.1 pc to ~1 pc) derived from the dichroic po- 
larization are running at a position angle of ^120° 
projected on the sky. 

4. Imaging polarimetry in the Ks band shows the dom- 
inant knotty nebulae of the shocked H2 emission is 
also polarized. They are due to dichroic absorption 
within the NGC 2071 core region. The derived mag- 
netic field direction projected on the sky (~115°) is 
consistent with that from the polarimetry of point- 
like sources. 

5. The magnetic field direction derived from the near- 
IR polarimetry is also consistent with that inferred 
from the previous 850 /im thermal polarimetry 

6. The magnetic fields in the NGC 2071 region of 0.1-1 
pc scale are running roughly perpendicular to the 
direction of the large scale outflow. We argue that 
the field strength is too weak to align the outflow in 
the large scale field direction via magnetic braking. 

This research was partly supported by MEXT, 
Grant-in-Aid for Scientific Research on Priority Areas, 
"Development of Extra-solar Planetary Science" , and by 
grants-in-aid from MEXT (Nos. 16340061, 17204012). 
We thank the referee, David Weintraub, for his helpful 
comments. 
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Table 1: Aperture photometry and 

polarimetry of IRS in the core region. 



IRS 


RA 


DEC 


J 


H 


Ks 


P(.') 


f)(J) 


P(H) 


0{H) 


p(Ks) 


<j{I< s ) 


Aperture 




(") 


(") 


(mag) 


(mag) 


(mag) 


(%) 


(°) 


{%) 


(°) 


(%) 


(°) 


(") 


IRS1 








16.51±0.02 


13.53±0.01 


11.31±0.01 


23.8±2.2 


178±4 


27.6±0.2 


177±3 


26.1±0.2 


177±3 


6 


IRS2 


8.7 


7.4 


14.54±0.01 


12.45±0.01 


10.52±0.01 


17.2±0.4 


155±3 


14.0±0.1 


159±3 


6.7±0.1 


155±3 


6 


IRS3 


-1.9 


5.6 


18.25±0.08 


15.51±0.01 


12. 19±0.01 






15.4±1.1 


128±4 


19.1±0.4 


124±3 


6 


IRS4 


5.0 


18.0 


15.03±0.01 


12.49±0.01 


10.86±0.01 


1.5±0.6 


37±10 


1.9±0.1 


115±3 


3.4±0.1 


121±3 


6 


IRS5 


3.7 


5.0 


16.83±0.02 


1 1.56±0.01 


12.39±0.01 


5.3±2.4 


141±12 


19.9±0.4 


115±3 


10.5±0.3 


150±3 




IRS6 


13.6 


27.3 


13.55±0.01 


11.37±0.01 


10.03±0.01 


4.0±0.2 


87±3 


3.6±0.1 


102±3 


3.6±0.1 


109±3 


6 


IRS7 


-12.4 


26.0 


18.70±0.12 


15.69±0.01 


12.44±0.01 






12.4±1.4 


17±4 


9.2±0.4 


20±3 


6 


IRSS 


-6.8 


5.0 


19.49±0.14 


16.62±0.01 


14.24±0.01 






9.5±1.7 


16±6 


7.7±1.2 


46±5 


3 


IRS8a 


-5.2 


4.3 


18.94±0.09 


16.65±0.01 


14.02±0.01 






6. Oil. 7 


16±9 


4.3±1.0 


163±7 


3 
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(b) Polarized Intensity (JHK 3 ) 
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Fig. 1. JHKs color composite images of intensity (a) and polarized intensity (b) of the NGC 2071 region. J is 
blue, H is green, and Ks is red. Field-of-view is 7 .'7 X 7f 7 in the left images, while it is 4' X 4' in the 
right images, (c) Notable IR sources and IR reflection nebulae (IRNe) are indicated. Field-of-view is 7.'7 X 7.'7. 
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5 h 47 m 7?0 6 S .5 6 5 .0 5 5 .5 5 s .O 4 S .5 4 s .O 3 5 .5 

RIGHT ASCENSION ( J2000 ) 

Fig. 2. KsHJ band polarization vector maps and Ks band identification map of the central 

N2071TR cluster region snnernosed on the intensitv images. Notable TR sources and TR 
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RIGHT ASCENSION (J2000) 

Fig. 4. Ks band polarization map of point-like sources superposed on the intensity image. Note 
that the polarization of the central source of IRN6 seen in the north-west, has a clear 
intrinsic polarization, which is almost perpendicular to the polarization of other point-like sources. 
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Fig. 5. KsHJ band histograms of position angle of point- 



;-likc sources. 
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Fig. 6. Comparison of the distribution of position angles at 1.63 (im and 850 fj,m. The 850 (im data are rotated by 90° 
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Fig. 7. Correlation between degrees of polarization in the H band and H — Ks color. The broken line shows a linear fit with a slope of 
6.0. 



